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Depolarization of quantum fields is handled through a master equation of the Lindblad type. The 
specific feature of the model is that it couples dispersively the field modes to a randomly distributed 
atomic bath. The depolarizing dynamics emerging from this approach is analyzed for relevant states. 
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I. INTRODUCTION 

Polarization of light is a fascinating concept that has 
deserved a lot of attention over the years. Today, the field 
is witnessing a renewed interest because of the fast devel- 
opments in optics, both on the applications and on the 
more fundamental aspects. Moreover, the polarization of 
a field is a robust characteristic, which is relatively sim- 
ple to manipulate without inducing more than marginal 
losses. Therefore, it is not surprising that many experi- 
ments at the forefront of quantum optics use polarization 
states [llllSli. 

The Stokes parameters provide, perhaps, the most con- 
venient description of polarization in classical optics and 
can also be extended in a natural way to the quantum 
domain |^ Q. They have the drawback in that they 
measure only second-order correlations of the field am- 
plitudes: while this may suffice for classical problems, 
higher-order correlations are important for the realm of 
quantum optics. This is the basic reason why the Stokes 
parameters do not distinguish between quantum states 
having quite different polarization properties 0, 1^ ■ 

Despite these subtleties, the Stokes parameters still 
provide the most commonly used definition of a degree 
of polarization |3, [lO, [Oi Ulli as we discuss in Sec. E] 
The degree of polarization is often desired to reach its 
maximum value, as well for classical as for quantum 
communication. The term depolarization has come to 
mean the decrease of the degree of polarization of a light 
beam when traversing an optical system. Our intuition 
strongly suggests that depolarization is due to decoher- 
ence. Roughly speaking, two basic mechanisms have been 
proposed to explain classically this phenomenon |l3| : 
decorrelation of the phases of the electric field vector and 
selective absorption of polarization states. The first cause 
is far more interesting to examine when the microscopic 
mechanisms responsible of depolarization are analyzed. 

To deal with this decoherence, some ingenious and 
practical tricks have been proposed, such as coupling 
each mode to a beam splitter where it is mixed with 
a vacuum [lj|. We wish to look at this problem from 
a more fundamental point of view. Of course, one can 



find in the literature a number of pure dephasing mod- 
els in which the environment is represented as a scat- 
tering process, which does not lead to a change of pop- 
ulations ^^ [13, [i3j HM H^' -'^o^ many applications, a 
good knowledge of dephasing is of utmost importance. 
This holds most prominently for quantum information 
processing, where the operation completely relies on the 
presence of coherence. It is hardly surprising that a 
lot of attention has been paid to dephasing in systems 
such as quantum dots [2J, |2l|, 123, Q , Josephsonjunc- 
tions |2J, |23j l2y| , or general quantum registers [23| , to 
cite only a few relevant examples. 

However, the study of depolarization is not only a ques- 
tion of pure dephasing dynamics. Polarization has an 
additional su(2) invariance that leads to a natural struc- 
ture of invariant subspaces. On physical grounds, we 
argue that this structure must be preserved in the evolu- 
tion, which makes previous models fail in this case. The 
main goal of this paper is precisely to provide a simple 
approach that goes around this drawback and provides a 
picture of the mechanisms involved in the depolarization 
processes. Since in classical optics the randomization is 
produced as light propagates in a material medium, our 
main idea is to couple the field modes dispersively to a 
randomly distributed atomic bath: the resulting master 
equation has a quite appealing structure that is exam- 
ined in some detail in Sec. lIIII for the case of single-mode 
fields. Finally, in Sec. IIVI we extend the formalism to 
multimode fields, while our conclusions are summarized 
in Sec. El 



II. POLARIZATION AND SU(2) INVARIANCE 

In the interest of retaining as much clarity as possible, 
we first recall some well-known facts about the polariza- 
tion structure of quantum fields. Since the formalism 
applies to fields of arbitrary wavefronts and frequencies, 
we consider a free transverse electromagnetic field with 
m spatiotemporal modes, whose positive-frequency part 



j=l s=± 



fkui 



E(+)(r,t) ^ iJ2J2 \h?r "-^^ Ujs(r)exp(-iwji) 



2eo 



(2.1) 

In the Dirac quantization scheme the modes are plane 
waves 



Uis(r) = -7^ ejs exp(ikj • r) 



(2.2) 



defined in some large volume V (which may be taken to 
be infinity later). For definiteness, we shall work in the 
base vectors e^g corresponding to circular polarizations 
(and s takes the values ±). The annihilation and creation 
operators ajs and a,^ obey the familiar commutation re- 
lations 



{ajs,a],J =Sjj,Ss 



(2.3) 



As pointed out by Karassiov |23| , there are specific ob- 
servables that characterize proper polarization properties 
of fields of the form H2.1(l . They correspond to the gener- 
ators of the group SU(2) of polarization gauge invariance. 
In the circular polarization basis they can be expressed 
as 






(2.4) 



They indeed furnish a Schwinger representation of the 
su(2) algebra: 



[Jz,J±]=±J±, [J+,J_]=2J, 

The Casimir operator is 

1, . , , . . N /N 



(2.5) 



J' - J: + i^iJ+J- + J-J+) = ^ hr + 1 ■ (2-6) 



The operator 






(2.7) 



represents the total number of photons and satisfies 

[iV,J]=0, (2.8) 

where J = (J^^, Jy, Jz), with J± = J^ i: iJy The total 
Hilbert space splits in this way in invariant subspaces of 
dimension A'^ + 1. Since there is no risk of confusion, 
we denote by the same letter the operator (|2.7|l and its 
eigenvalue, which is the total number of photons. 

The quantities J are then measurable in photon- 
counting experiments |29j. In fact, for a single-mode field 



they coincide up to a factor 1/2 with the Stokes opera- 
tors pMlj. It is thus natural to extend this identification 
and define 



-^« 



(2.9) 



for TTi-mode fields. This allows one to parallel the classical 
idea and define the degree of polarization as 



P 



{N) 



(2.10) 



Note that s — {S)/{N) is the polarization vector in clas- 
sical optics, which defines the Poincare sphere. Equa- 
tion (|2.1U|) can then be interpreted as the distance from 
the point represented by s to the origin, which is associ- 
ated with the unpolarized light. 

To highlight the SU(2) invariance of p olarization, it 
has been proposed by several authors [sj, Hj, Hj, Hj, 113 
to define unpolarized light as the field states that remain 
invariant under any linear polarization transformation, 
which in experimental terms can be accomplished with 
a combination of phase plates and rotators (that pro- 
duce rotations of the electric field components around 
the propagation axis). Any state satisfying the invari- 
ance condition will also fulfill the classical definition of 
an unpolarized state, but the converse is not true. 

The density operator of such quantum unpolarized 
states can be written as 



gunpol = ^ rN Iat 



(2.11) 



N 



where Ijv is the identity operator in the invariant sub- 
space with A^ photons and all r^ are real and positive 
constants. The condition of unit trace imposes 



^(A^+l)rjv = l. 



(2.12) 



N 



The vacuum state is the only pure quantum state that is 
unpolarized, and the unpolarized mixed states are totally 
mixed in each subspace. Note that in each invariant sub- 
space, gunpoi is fully random: we claim that in any phys- 
ical depolarization process this fundamental structure of 
invariant subspaces must be preserved. 



III. LIGHT DEPOLARIZATION: SINGLE-MODE 
FIELDS 

A. Master equation for pure dephasing processes 

We begin by focusing our attention on a single-mode 
field coupled to a bath system. The total Hamiltonian is 



H = H, 



field 



H, 



bath 



H. 



(3.1) 



Here -ffficid is the Hamiltonian of the mode under consid- 
eration 



ffficld = ^ X! "«"^ ' 
s=± 



(3.2) 



while i?bath describes the free evolution of the environ- 
ment. We do not make any hypothesis on the precise 
kind of bath, and only assume that it is so large that its 
statistical properties are unaffected by the coupling with 
the system. 

In the interaction picture and under the usual weak- 
coupling and Markov approximations, the master equa- 
tion takes the general form 17] 



m 



1 



dT TrB{[i/int(i), [Hintit-T), g{t)(g>gB]]} , 

(3.3) 

where Tre is the partial trace over the bath variables and 
£•101(0 ^ £'(0 ® £'B is the density operator for the total 
system. In addition, we take the mode and the bath 
initially independent. Henceforth, g{t) will denote the 
reduced density operator for the field. 

A common way to couple the system to the bath is 
through the basic interaction Hamiltonian 



-ffint — fi 



A s=± 



{k\ Txal + k\ V\as) , 



(3.4) 



where nx are coupling constants, T\ and V\ are annihila- 
tion and creation operators for bath quanta and the sum 
over A runs over all the accessible bath modes. Note that 
a photon can lose energy by creating a bath quantum, 
and conversely. 

We assume for simplicity that the bath is at zero tem- 
perature, which is tantamount to neglecting stimulated 
processes. The master equation (|3.3|) is then 



^ = E?^ 



(3.5) 



where L[Cs\ are the Lindblad superoperators [3a| 

C[Cs] g - 2CsgCt ~ {C]C„ g} . (3.6) 

and "fs are the decoherence rates associated with each 
polarization. As was recognized earlier ^J^ ,3SJ , this form 
of the master equation is completely positive, which is 
the only dynamical evolution of an open system ensuring 
that the state of the system remains physically valid at 
all times. 

Equation 13.5|l corresponds to two standard indepen- 
dent Lindblad decaying processes, one for each basic po- 
larization. Since the coupling l|3.4|l transfers energy be- 
tween the system and the bath, the net result is a damp- 
ing of the beam intensity. In fact, one can check that, 
irrespective of the initial conditions, the stationary state 
is always g{t ^ 00) = |0,0>(0,0|, where |0,0) = |0) + |0)„ 
is the vacuum state for both circular polarizations. Ob- 
viously, this cannot describe depolarizing processes. 

As we have anticipated in the Introduction, depolar- 
ization implies a pure dephasing process. In mathemati- 
cal terms, this means that the interaction must commute 
with the system Hamiltonian, so no energy is transferred 
and only the phase changes. Various theoretical scenar- 
ios have been proposed to that end. Apart from minor 



details, all of them can be modeled by an interaction 
Hamiltonian of the type [l7| 

-ffint = X! X! (^-^ ^\a,las + kI T{alas) , (3.7) 

A s=± 

which can be viewed as a scattering process in which 
a bath quantum can be absorbed or emitted, but the 
number of photons in each polarization is preserved. 

Assuming again a zero-temperature bath, the master 
equation for (|3.7|) can be written as 



g=Y^ -^C[alas] g. 



(3.8) 



To illustrate the physics embodied in this equation, let 
us focus on the fully quantum regime and consider a one- 
photon state. We are then within a two-dimensional in- 
variant subspace, whose basis will be labeled 



|+> = |1> + |0)_, |-) = |0) + |1)_ 



(3.9) 



The density matrix in this two-dimensional subspace can 
be expressed as 



W 1 + s. s^ *^y 1 =_(£ + s-cr), (3.10) 



2\s 



~x ~r t'Jy 



tSti 1 Sz 



where s is the polarization vector defined after Eq. (|2.1(J|) 
and can be written as s = 2Tr(£»cr), cr being the Pauli 
matrices. In this subspace the model (|3.8|l can be easily 
recast in terms of s; the final solution reads as 



sy{t) = s,(0)e-(^++^-)*/2, 
Sz{t) - s,(0). 



(3.11) 



This interaction indeed preserves the invariant subspaces, 
but produces no thermalization: its effect is merely to 
maintain the occupation probabilities, while erasing all 
coherences. Although this can model appropriately many 
interesting phenomena, it fails to describe depolarizing 
effects, because depolarization leads to a uniform dis- 
tribution of photons in each invariant subspace and not 
only to a pure dephasing. In other words, depolariza- 
tion not only preserves the invariant subspaces, but the 
steady state in each one of then must be a diagonal state. 
We believe that these conditions are essential to ensure 
a correct description of depolarization. 



B. Modeling depolarization with a nonresonant 
randomly-distributed atomic medium 

Apart from the previous drawbacks, generic models of 
the type H3.7|l suffer from the problem that one cannot 
provide a physically feasible mechanism leading to an in- 
teraction of that form. To solve these difficulties we shall 
take another route. We recall that in classical optics, the 



electric field of the light experiences a decorrelation as it 
propagates through a material medium. If a± denote the 
classical field-mode amplitudes, such a decorrelation can 
be taken into account by assuming that a± — *■ a±e"^'^ , 
where (p± are random phases. 

We wish to explore this picture in the quantum world. 
In consequence, we assume that the field propagates 
through a material medium represented by a collection 
of two-level atoms. Since there is no net energy transfer 
between the field and the atoms, their interaction must 
be necessarily dispersive. In short, our basic system is 



Hgys — -ff field + -ffat + V , 



(3.12) 



where 



^at = '^-hujxcrl, 



(3.13) 



V = h 






idXs 0-A < + g'L f^A Os) 



The form of the interaction V assumes, as it is usually 
done, that the atoms interact with the field mode in the 
dipolar and rotating-wave approximations. In addition, 
the atoms are randomly distributed so the coupling con- 
stants gxs have random phases 39] . 

If A\ — LU\—uj is the detuning, we must consider the far 
off- resonant regime \g\s\ ^ A^. In such a limit, we can 
adiabatically eliminate the nonresonant interactions in 
Eq. (|3.13|) and obtain the following effective Hamiltonian 
that describes our system ^] : 






1 



A 



A, A' S — ± 



E E \^9xs9ys (^^ + ^ j ^t<^x' . 

(3.14) 



where 

M=Y, alas + Y, -ai 

.s=± A ^ 

is the conserved excitation number operator and 

Jx 



(3.15) 



r!> 



Aa + ^A + 



Aa' 



(3.16) 



Jx = 2\gx\^{N + J+e''^>- + J^e-''^-). 



Here we have written gx± = IsaIc^*'''^^^, with (px being 
random phases. This refiects the fact that the relative 
phase between the atomic dipole and the field mode is 
random. The parameter 



^A = 



\9x\ 



(3.17) 



represents just a small frequency shift. 

Moreover, it is well known that the atoms decay irre- 
versibly. This usually assigned to their interaction with 
the continuum of modes of an additional thermal electro- 
magnetic environment. In such a case, the density matrix 
for the system H3.12II evolves according to J17l | 



£'sys(0 



l.-'^sysi ^sys 



A 

(3.18) 



where 7a is the decay constant of the Ath atom due to 
its coupling to the thermal environment with fix excita- 
tions. As is well established J4l| 142, |4J| , the properties 
of a random medium are well reproduced in the high- 
temperature limit: fi-A ^ 1. In this limit, the effect of 
spontaneous emission can be disregarded in comparison 
with the stimulated emission processes. Emission into 
the reservoir and absorption from the reservoir therefore 
become identical; i.e., they balance each other in the sta- 
tionary state: the emission and absorption processes de- 
pend solely on the initial population of the atomic state. 
Consequently, the steady-state reduced density operator 
is approximately given by a mixture of equally popu- 
lated atomic states. The density matrix of the Ath atom 
thus becomes diagonal (^ia — ^^) ^'^d the effect of the 
last term in Eq. (|3.14|) is negligible and we shall omit 
it henceforth. The effective Hamiltonian for our system 
then reads as 



H. 



hujN- 



Y.\ 



Mlxu 



Af^A- 



(3.19) 



As indicated in the Appendix, in this limit we can also 
adiabatically eliminate the atomic variables and obtain 
a master equation that, after averaging over the random 
phases, reads as 

e^~i[LoN,B] + lc[N]B+^-C[J+]g+^C[JAQ, (3.20) 

where g(t) = Trat[psys(i)] is again the reduced density 
operator for the field mode and the decoherence rate 7 is 



7 



4E 



I.9AI 



7AA^nA 



(3.21) 



Equation (|3.20() is our central result. We observe that it 
preserves the SU(2) invariant subspaces and the steady 
state in each iV-photon subspaces is a completely random 

state 



Q{t ~* 00) = 



1 



N+l' 



■ a. 



(3.22) 



This can be considered as the major advantage of our 
approach. 

A rough estimate immediately shows that the depo- 
larization rates (|3.21|) appearing in this model are very 



small when compared with other typical system parame- 
ters, which is in agreement with the experimental obser- 
vations. 

The terms /^[Ji] g describe depolarization in each in- 
variant subspace, meanwhile the action of £[iV] g therein 
is trivial. Nevertheless, this C[N] g is responsible for the 
relative phase decay between blocks of the density matrix 
corresponding to different excitation numbers. 

To compare with the previous discussion in Eq. H3.11|l . 
we consider again the one-photon case. By recasting the 
depolarization master equation l|3.2U|l in terms of s, one 
obtains the solution as 



s4t) = s.(0)e~^*, 
sy{t) = Sy{0)e-''\ 



(3.23) 



which displays the aforementioned desirable properties. 
The degree of polarization (|2.1(J|I of this mode evolves 
then as 

P{t) = [4(0) + 4(0) + 52(0)6-2-^^/^6-^* . (3.24) 



IV. LIGHT DEPOLARIZATION: MULTIMODE 
FIELDS 

For the multimode case, it seems reasonable to assume 
that the depolarization is independent in each mode. The 
corresponding master equation can be thus simply writ- 
ten as 



Q = -^^^][Nj^Q\ 



i=i 



+ Eyi^t^j-] s + '^[Jj+\ e + '^[Ji-] rf- (4-1) 

One can check that the properties discussed in the pre- 
vious section also hold in this multimode case. 

For concreteness, henceforth we restrict our attention 
to the case of a two-mode field. We have then a repre- 
sentation of SU(2) (X) SU(2), and the whole Hilbert space 
can be decomposed in irreducible subspaces. For the case 
where there is one photon in each mode, we have, using 
the standard terminology of SU(2) representations, that 
D^/^ ^ jji/^ = £)0 £)i^ where D^ is the invariant sub- 
space with eigenvalue j. 

The one-dimensional subspace D^ is spanned by the 
singlet state 

IV'-)-^(l+)i|-)2-|-)i|+)2), (4.2) 

while D^ is spanned by the triplet states 

IV'^) = ^(I+>1 |->2 + |->1 |+>2), |+>1 |+>2, |-)l |->2 . 

(4.3) 



Here the states |±) for each mode are given as in H3.9|l and 
the subscripts 1 and 2 label the corresponding modes. 

We take our system to be initially in a generic density 
matrix such as 



P(0) 



gii(O) ei2(0) ei3(0) giM 

e2l(0) g22(0) ^23(0) P24(0) 

e3l(0) ^32(0) ^33(0) ^34(0) 

e4l(0) ^42(0) ^43(0) ^44(0) 



(4.4) 



(4.5) 



where we have employed the standard eigenbasis 

|l) = |+)l|+)2, |2) = |+)ih)2, 

|3) = |-)l|+)2, |4) = |-)lh)2. 

After some computation, the solution of the master equa- 
tion ()4.1|l can be shown to be 

g2i{t) - i[g2i(0)(l + e-2'^^*) 

+ £<43(0)(l-e-2'n*)]e-72t^ 

gziit) = i[g3i(0)(l + 6-2''^*) 

+ e42(0)(l-e-272t)]e-7l*^ 

gii{t) = e4i(0)6-(^i+^^)*, 

g^2{t) - e32(0)6-(^^+^^)*, 

^42(0 = ^[g42(0)(l + e-2''^*) 

+ e3i(0)(l-e-272*)]e-7i*^ 

^43(0 = ^[to(0)(l + 6-271*) 

+ e2i(0)(l-6-2''i*)]e-7^*, (4.6) 

for the nondiagonal elements, and 

giiit) = i(l + {2[f.„(0) + e22(0)]-l}e-27i* 
+ {2[pn(0) + e33(0)]-l}6-272* 

+ {2[gn(0) + e44(0)]-l}6-2(7i+7.)*)^ 

te(i) = i(l + {2^2(0) + eii(0)]-l}e-27it 

+ {2[e22(0)-fe44(0)]-l}6-272t 

+ {2[te(0) + e33(0)]-l}6-2(7^+^^)*), 

^33(0 = i(l + {2^3(0) + e44(0)]-l}e-27i* 

+ {2^33(0) -)-eil(0)]-l}6-27=* 

+ {2^3(0) + ^22(0)] - l}6-2(7i+7^)*) , (4.7) 

for the diagonal ones. Here 71 and 72 are the decoherence 
rates for modes 1 and 2, respectively. Note that the trace 

condition gives ^44 = 1 — pn — £122 — e33- 

In this basis, the operators J defined in Eq. H2.4(l . take 
the explicit form 

Q -i -i 

-'y- 2 \ i -i 












(4.8) 



J. 



10 


0-1 



In consequence, the degree of polarization 12.10|l is now 

P{t) = [|ei2(i)+£'l3(i)+e24W + g34(t)|' 

+ \Qii{i)-QAA{t)\^]. (4.9) 

This is a quite compact resuh. Consider, for example, a 
pure disentangled state such as any element of the basis 
(14.5(1 . say q{Q) = |1)(1|. Using H4.6|l and ((4.7|l we imme- 
diately get 



P{t) 



{e-^'^^' + e-^''^*) 



(4.10) 



On the other hand, for the maximally entangled Bell-like 
states 14411 



|^±)(^^ 



we have 




1(0 1 ±1 
2 I ±1 10 





P{t) = -e-(''i+''2)* . 



(4.11) 



(4.12) 



Now, owing to the strong correlations, the decoherence 
rates for both modes always appear together. We stress 
that even the singlet state will become depolarized. This 
is due to the fact that the depolarization mixes differ- 
ent invariant subspaces. Note also that although these 
two examples are pure states, they show quite a different 
evolution of the degree of polarization. 



APPENDIX A: DERIVATIO N OF THE MASTER 
EQUATION (IST^ni) 

Let us define the following operators 



Aig^QxQ: 



A^g = g^x , 



(Al) 



which satisfy the following commutation relations 



[iV±,M,] = TN^Sxy 
[iV±,7Vj:,] = TNtSxx' 
m.Af] ^ 0, 



(A2) 



where a,/? = ±,Z,r. Apart from their interesting alge- 
braic properties, about which we are not concerned in 
this paper, the operators I(A1(I allow one to recast the 
master equation (|3.18l) in a more appropriate form. In 
fact in the high-temperature limit we can replace -ffsys by 
the effective Hamiltonian ((3.19|l and we have 

&ys = -y^7A»Agsys 
A 

+ J2^jxnxiN+ +N-)^ 2zA{Nl + 2zAlNl]g,y, . 

A 

(A3) 

The dissipative term in this equation is much bigger 
than the Hamiltonian term. This allows us to adiabat- 
ically eliminate the atomic degrees of freedom. We can 
e.g. make use of the perturbative method proposed in 
Refs. 45] and 46] (see Ref. !47i| for a complete account 
of the application of the method to effective master equa- 
tions) . This requires us to represent Eq. I(A3(I in terms of 
diagonal and raising-lowering operators. To this end, we 
first diagonalize the dissipative part of l(A3|l by applying 
the following 7r/2 rotation 



V. CONCLUDING REMARKS 

What we expect to have accomplished in this paper is 
to present a comprehensive theory of the depolarization 
of quantum fields. In our model the field modes couple 
dispersively to a randomly distributed atomic bath: the 
resulting master equation has unique properties that we 
have explored in detail. 

We hope that the tools introduced here could be of 
interest in treating depolarization in the fully quantum 
regime, especially for biphotons, which have strong impli- 
cations in areas of futuristic technologies such as quan- 
tum computing, quantum cryptography, and quantum 
communications. 



V 



exp 



2^2^- 
A 



iv: 



(A4) 



which leads to 

^-^y- = -5IbA^A+7A"A(iVI + M)/2 

A 

- i{A\~A\){Nt+N^) 
+ i{A\ + A\){N{-Nl)]g^y^. 

Next, we apply the following "small rotation" 



U = exp 



, , AT 



AK 



'^ixnx 



m 



n: 



(A5) 



(A6) 



After a lengthy but otherwise straightforward calculation 
and applying the inverse of the transformation (|A4() one 
finally gets 



Psys 



-Ei^^"^[l-^A+-^A" 



A 

»(^A 






N{)}Qsys 



A 



27A?^A 



A{ 



27A?^A 



(iV++7V-)e,y,(A7) 



where the field mode is expressed in a rotating frame. 
This transformed equation has the virtue of containing 
only diagonal terms so we can immediately trace over 
atomic variables to obtain the master equation for the 
density matrix of the mode g(t) = Trat[esys(i)]- Since 
each atom in the medium is in a statistical mixture, gx — 
il, the first term in (|A7|) vanishes and one obtains 



E 



1 



27AAAnA 



^Jx]q, 



(A8) 



where J'\ has been defined in Eq. (|3.16() . If we average 
over all the random phases ipx we get 



g^lc[N]g+^C[J+]g+^C[J-]g, (A9) 



with 7 given in (|3.21|) . This is precisely the master equa- 
tion (|3.20() . written in the interaction picture. 
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